Reactor Antineutrino Anomaly with known #13 
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We revisit the reactor antineutrino anomaly using the recent reactor flux independent determina- 
tion of sizable #13 by considering the full set of the absolute reactor v K flux measurements. When 
normalized to the predicted flux of Mueller et al. pj , the new world average, after including results 
from Palo Verde, Chooz, and Double Chooz, is 0.959 ± 0.009 (experiment uncertainty) ± 0.027 (flux 
systematics) . Including the data with kilometer baseline, the new world average is only about 1.4a 
lower than the unity, weakening the significance of the reactor antineutrino anomaly. The upcoming 
results from Daya Bay, RENO, and the Double Chooz will provide further information about this 



I. INTRODUCTION 

The term "reactor anomaly" was coined by Mention et 
al. [H who noted that the average of the experimentally 
determined reactor antineutrino flux at reactor-detector 
distances < 100 m accounts for only 0.943 ± 0.023 of the 
reevaluated theoretical expectation of Ref. [l[ . In addi- 
tion to the 19 experimental results obtained with detec- 
tors distant less than 100 m from the reactor source, we 
are able now to include in the analysis also the results of 
the Chooz SB and Palo Verde @ as well as of the Dou- 
ble Chooz [6|— lS|] experiments, where the detectors were 
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further away from the reactor complex. In these cases, 
the corresponding experimental results need be corrected 
for the flux loss associated with the known value of the 
mixing angle #13, which was determined in a model in- 
dependent way by comparing the count rates in two es- 
sentially identical, but separated in distance, detectors. 
We use the value sin 2 28 13 = 0.089 ± 0.011 obtained in 
the Daya Bay experiment and confirmed by the 

RENO experiment fl2l ]. The corresponding correction is 
easy to apply using the formula for the survival proba- 
bility 



P, 



1 - sin 2 2^i3 (cos 2 6>i2 sin 2 A31 + sin 2 6>i 2 sin 2 A 32 ) 



6»ia sin 2 26> 12 sin 2 A 21 , 



(1) 



with A, 



'•3 — I ~~ l-27|Amfj| E(MeV) ■ 

ing angles and mass-squared differences used in the sim- 
ulation are taken from Ref. fl3j assuming normal mass 
hierarchy (Am^ = Am| 2 + Am^): 



L(m) 



Values of mix- 



sin 2 20i2 = 0.857 ± 0.024 

Am 2 : = (7.50 ±0.20) x lO^eV 2 
Am 2 32 = (2.32 ± 0.12) x 10~ 3 ey 2 



II. ANALYSIS DESCRIPTION 



(2) 



In the following, we will explain in detail the inputs to 
our analysis as well as the chi-square method. For the 
19 experimental results obtained with detectors distant 
less than 100 m from the reactor source, the measured 
fluxes normalized to the prediction of Ref. [l[ (noted as 
"ratio") are taken from Ref. Q and tabulated in Table. Q] 



I 

together with the detector technology, fission fractions, 
distance to the reactor core, and year of publication. The 
average survival probability P SU r for each experiment is 
calculated by integrating over the neutrino antineutrino 
spectrum [l| convoluted with the inverse beta decay cross 
section with sin 2 29i 3 = 0.089. There are two uncertain- 
ties listed in Table HI The a err represents the total un- 
certainty on the ratio. The a corr represents the part of 
uncertainty correlated among different experiments. In 
particular, there is a 2.7% uncertainty in a corr coming 
from the uncertainty in the predicted reactor flux. (In 
first versions of Ref. [2J , the reactor flux uncertainty was 
assumed to be 2.7%, which was replaced by 2% in the 
final version. In this analysis, we chose the more conser- 
vative number (2.7%).) Since it will not affect the relative 
differences among different measurements, we define the 
reduced uncertainties: 
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tainty. We change the <r c „, , 

II [3 from 3.7% (originally quoted in Ref. 0) to 2.7% 
(the reactor flux uncertainty only). The original 3.7% 
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Experiment # 

FIG. 1. (color online) The covariance matrix of the reduced 
experimental uncertainties. The number of experiments can 
be found in Table. [T] The off-diagonal terms show the corre- 
lation among different experiments. 



assumes almost full correlation between SRP-I and SRP- 
II, which can not explain the apparent differences be- 
tween the two ratios (0.952 vs. 1.018), indicating less 
correlation between the two experiments. For the same 
reason, we reduce the a corr between the ROVN088-1I 
and ROVN088-2I from 6.9% to 5.7%, which is the final 
reported uncertainty for ROVN088 experiments jiEj . 

There are three reactor cores in the Palo Verde ex- 
periment The distances between detector and each 



Here, W 1 is the inverted covariance matrix. The vector 
R contains the reported ratios from all 23 experiments 
(tabulated in Table. Q} . The absolute normalization ratio 
r is treated as a free parameter. The vector P sur contains 
the predicted average survival probabilities given a value 
of sin 2 26*13. The values of P SU r using sin 2 26*13 = 0.089 



To be consistent with other experiments, we choose the rate- 
only sin 2 2$i3 results, which represent simple measures of the 
disappearance in the total number of events. 



reactor core are 750 m, 890 m, and 890 m. In calcu- 
lating the average survival probability P sur , we assume 
that all three reactor cores have equal power. The result 
is compared with Pj™ m assuming full power in only the 
750 m reactor and Pf^ m assuming full power in the 890 
m reactors. The differences are quoted as an additional 
uncertainty, which is only about 5% of the total reduced 
experimental uncertainty. 

There are two reactor cores in Chooz experiment 0, HJ . 
The distances between the detector and each reactor core 
are 998 m and 1115 m. Similar procedure is applied to 
calculate the uncertainty for the equal power assumption. 
The resulting uncertainty is about 6.2% of the total re- 
duced experimental uncertainty. The fission fractions are 
assumed to be the same as those from Double Chooz [3] ■ 
We also calculated the average P sur by varying these fis- 
sion fractions. The differences are negligible. 

The Double Chooz experiment is conducted at the 
same location as Chooz. With a single detector, the re- 
cent rate-only analyses of the data from delayed neutron 
capture on Gadolinium (n-Gd) and delayed neutron cap- 
ture on hydrogen (n-H) reported the value of sin 2 26>i3 = 
0.170 ± 0.052 and sin 2 2(9i 3 = 0.044 ± 0.060 1 0, by 
anchoring to the short-baseline Bugey-4 results [lg|, re- 
spectively. Although the measured flux normalized to 
the prediction of Ref. [l[ has not been reported, we can 
deduce such ratios using the reported fission fractions 0] , 
the reported values of sin 2 29i 3 @, [1] , and the Bugey-4 
results. The a r J.f uced is dominated by the uncertainties of 
reported sin 2 26*13, with additional uncertainties coming 
from the equal power assumption. The a r c ^ ced are calcu- 
lated from the reduced experimental uncertainty a^ duced 
from Bugey-4. Furthermore, there arc additional corre- 
lated uncertainties between the n-H and n-Gd measure- 
ments due to the equal power assumption. The final co- 
variance matrix W using reduced uncertainties is shown 
in Fig. [U 

The x 2 function used in this analysis is constructed as 
follows: 



(4) 

I 

are tabulated in Table. Q] The last term in Eq. (j4|) repre- 
sents the constrain on sin 2 26*i3 from the latest Daya Bay 
results [HI- 

III. RESULTS 

In the following three figures, we show the results of 
all 23 measurements and the deduced ratios after mini- 
mizing the x 2 defined in Eq. The global average is 
determined to be 0.959±0.009. In FigJU the results are 
shown in analogous way as in Ref. [2|], i.e. as a func- 



X \r,sm' 20i 3 ) = (r • P sur (sin* 29 13 ) - R) T W-\r ■ P s „ r (sin 2 29 



R) 



(sin 2 26*13 - 0.089) 2 
0.011 2 
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tion of the distance from the corresponding reactor core. 
We combine results at same baseline together for clarity. 
The corresponding x 2 /Dof = 23.8/22. We stress that 
our error bars do not include the reactor flux uncertainty 
(2.7%), hence they appear smaller than those in 0]. 



Distance 



o 1.2 
+■» 

(0 

r 

1.1 
1 

0.9 
0.8 
0.7 
0.6 



TT 



10 



10 2 



— — World data 
World Average 

I 1-a Exp. Unc. 

I 1-a Flux Unc. 

, i 

10 3 



Distance (m 



FIG. 2. (color online) The reactor v e capture rate as a 
function of the distance from the reactor, normalized to the 
theoretical flux of Ref . [l[ ■ The horizontal bar represents the 
global average and its la error bar. The 2.7% reactor flux 
uncertainty is shown as a band around unity. We combine 
results at same baseline (e.g. Chooz, Double Chooz n-H and 
n-Gd results) together for clarity. 
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FIG. 3. (color online) The same 23 experimental results 
as in previous figures are plotted by the detection method 
employed. See captions of Fig. [2] for details. 



related uncertainty between SRP-I and SRP-11 from 3.7% 
to 2.7% (the reactor flux uncertainty only), since the orig- 
inal 3.7% assumes almost full correlation between SRP-I 
and SRP-II, which can not explain the apparent differ- 
ences between the two ratios (0.952 vs. 1.018). With 
fixed total experimental uncertainties, this change effec- 
tively increases the significance of SRP experiments and 
leads to about 1% larger world average. 

In addition, our results are larger than those reported 
in Ref. [ljj and Ref. [l8j], which also include the kilo- 
meter experiments with known #13. The result reported 
in Ref. pj} included the Gallium neutrino data fl9j - [25lj . 
which was not included in our reactor antineutrino anal- 
ysis. They also did an analysis by including a RENO pre- 
liminary result from the absolute flux analysis. However, 
such analysis of the RENO experiment has not been, to 
our knowledge, released and is not finished as yet [26| . 
In Ref. [l8j], the measured experimental fluxes are nor- 
malized to the predicted flux of Huber [13] with a new 
neutron lifetime 881.5s (2011 update of PDG [Hj]). The 
change in reactor flux model and the neutron lifetime 
leads to in average 1.6% lower ratios than what we used 
in this work (tabulated in Table. HI . The rest of differences 
come from the treatment in the correlated uncertainty of 
SRP experiments and the uncertainty of the reactor flux 
prediction (2% used in Ref. 18] vs. 2.7% used in this 
work) . 

One of the main purposes of this work is to illus- 
trate the impact of kilometer experiments to the results 
of Mention et al. @. Therefore, we have adapted the 
same neutron life time (885.7s) as in Mention et al. 0. 
The current recommended neutron lifetime from the 2012 
Particle Data Group [l3[ is 880.1s. Using the latest neu- 
tron life time would lead to about 0.63% reduction in the 
average ratio. 
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The new global average is somewhat larger than the 
0.943 value of Ref.Q quoted earlier, weakening the sig- 
nificance of the reactor antineutrino anomaly. There are 
two reasons for this difference. First, we include more re- 
cent and more distant experiments, of which Palo Verde 
and Chooz have larger rates. Second, we change the cor- 



FIG. 4. (color online) The reactor i> e capture rate as a 
function of the year when the measurement was published, 
normalized to the theoretical flux of Ref. [l[ . We combine re- 
sults at same year (e.g. Double Chooz n-H and n-Gd results) 
together for clarity. See captions of Fig. [2] for details. 
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FIG. 5. (color online) Best fits are shown with and without 
sin 2 2#i3 constrains from Daya Bay. The corresponding 68% 
(Ax 2 < 2.3) and 90% (Ax 2 < 4.6) confidence intervals are 
shown as well. 

Fig. [3] shows results from all 23 experiments again ar- 
ranged by the detector technology. Five different tech- 
nologies were used in the 23 experiments to record the 
v e capture on protons. In the ILL [28j and Goesgen ex- 
periments [29| the liquid scintillator targets cells were in- 
terspaced with the 3 He neutron counters. In the Bugey- 

3 experiments [3(1 [HJ the liquid scintillator was loaded 
with 6 Li to detect the neutron captures. In the Bugey- 

4 and Rovno91 (HJ only the neutron captures were 
detected. The detector in these two experiments con- 
sisted of a water target with embedded 3 He detectors. 
In the Krasnoyarsk |33j and Rovno88 flol ] experiments 
again only the total capture rates were measured. The 
detectors consisted of the polyethylene neutron moder- 
ator with 3 He neutron counters embedded in them. Fi- 
nally, the Savannah River experiments [lij , the Rovno88 
[loT ] and the Chooz, Palo Verde and Double Chooz exper- 
iments El use the Gd loaded liquid scintillators. The 
Daya Bay and RENO experiments are employing that 
technology as well. The results from the Gd loaded scin- 
tillator experiments are, with the exception of Double 
Chooz, in general higher than the new world average. 

We should also note that the experiments with Gd 
loaded scintillator were in general carried out at a later 
time. In Fig. |4]we show the same data, but now arranged 
as a function of the year when the corresponding mea- 
surements were published. One can see that the more 
recent experiments, except the Double Chooz, appear to 
have higher rates than the earlier ones. If such tendency 
is true, it could either be due to difference in detector 
technologies or improvement in data acquisition or anal- 
ysis methods. However, the x 2 of the global fit, x 2 /Dof 



2 In Double Chooz n-H analysis [j|], the neutrino events are largely 
generated in the liquid scintillator region without the Gd loading. 



= 23.8/22, signifies that the 23 reactor flux determina- 
tions are mutually consistent. Therefore, our observation 
might be simply due to statistical fluctuations. 

In addition, we also use these 23 experiments to ex- 
tract sin 2 20i3 by minimizing the x 2 in Eq. with- 
out the penalty term ( ^'" ^on" )• The best-fit 
sin 2 2#i3 is determined to be 0.084 ± 0.030 (one dimen- 
sion: Ax 2 =X 2 -xLmrrmm < !)■ The non-zero sin 2 20 13 
is dominated by the latest Double Chooz results, and is 
also consistent with the latest Daya Bay results [11] . 
Fig. [5] shows the best-fit sin 2 2£>i3 and ratio r with and 
without the Daya Bay constrains. The corresponding 
68% (two dimensions: Ax 2 < 2.3) and 90% (two dimen- 
sions: Ax 2 < 4.6) confidence intervals are shown as well. 
The correlation between the ratio and sin 2 29\z is rather 
weak due to the strong correlation between results from 
Double Chooz, which dominate the non-zero sin 2 2#i3 
extraction from these 23 experiments, and result from 
Bugey-4, which is the most precise short baseline mea- 
surement. Therefore, it is not surprising that the precise 
sin 2 2f5i3 value from Daya Bay, which is very close to the 
best fit sin 2 26>i3 value from these 23 experiments, does 
not improve the knowledge on the ratio. 



IV. SUMMARY AND DISCUSSION 

With the known 6*13, we include results from Palo 
Verde, Chooz, and Double Chooz, and re-evaluate the 
reactor antineutrino anomaly. The new world average 
is determined to be 0.959 ± 0.009 (experiment uncer- 
tainty) ± 0.027 (flux systematics), when normalized to 
the predicted flux of Mueller et al. [l|. The new world 
average is about 1.4c lower than unity, weakening the 
significance of the reactor antineutrino anomaly. We 
also show that the analysis of all 23 experiments yields 
sin 2 20i3 = 0.084 ± 0.030 in agreement with its indepen- 
dently obtained value from the Daya Bay and RENO 
experiments. 

The forthcoming absolute flux analysis of the Daya 
Bay, RENO, and Double Chooz experiments will clearly 
shed more light on this issue. The final answer to the 
question of agreement between the predicted reactor v e 
flux and the measured rate clearly depends on the results 
of these experiments as well as on the careful analysis of 
the predicted flux and its uncertainties. 

The most popular explanation of the anomaly, if it 
turns out that it is statistically significant, is the ex- 
istence of additional sterile neutrinos. The consistency 
of the rate of the short baseline neutrino experiments, 
demonstrated here and also e.g. in Ref. 0, [l8[ suggests 
that such additional neutrinos must have large enough 
masses that the corresponding oscillation length is at 
most few meters. Experiments sensitive to such short os- 
cillation length would be able to convincingly prove the 
existence of the sterile neutrino and allow determination 
of their masses and mixing angles. 



5 



# 


result 


Det. type 


235 t t 

u 


239p„ 
.r U 


238 TT 
U 


.r U 


ratio 


O err \ So ) 


&corr y/O) 




-L sur 


Year 


1 


Bugey-4 


3 TT.-, i TT ( \ 

rie + H2O 


0.538 


0.328 


0.078 


0.056 


0.942 


3.0 


3.0 


15 


U. 999987 


1994 


o 
z 




lie ri2 y) 


a fti a 
U.014 


a o7/i 
U.z / 4 


A (Y7 A 

U.U / 4 


A HQS 

U.Uoo 


A O/IA 

u.y4U 


o.y 


Q A 
O.U 


1 Q 

10 


A OOOOS1 

u.yyyyoi 


1 001 
ryyi 




Double Chooz 


(jCl-LiO 


a /iQft. 


U.oOi 


U.Uo / 


U.UOO 


A QftA 

U.oOU 


0. / 


Q A 
O.U 


OOQ 111C 

yyo-iiio 


A QKA 

U.yo4 


oAi 

zUlz 


OQ 

16 


Double Chooz 


Lb [n-rij 


A /I Aft 


U.ool 


A AC7 


U.UOO 


A AOA 

U.yzU 


A A 
4.U 


Q A 
O.U 


nnc 111c: 
yyo-rr ro 


A ClKA 

U.yo4 


oAi 
zUrz 


o 

o 


Bugey-3-1 


6t : t c 
Ll - Lb 


A COO 

U.000 


A QOO 

0.328 


A n7o 
U.U/8 


U.Uoo 


A A /I ft 

u.y4o 


A O 




lo 


A AAAAQT 


1 AAC 

iyyo 


/] 
4 


Bugey-3-II 


67 ■ T Q 


U.000 


U.ozo 


A A7S 

U.U / 


a 

U.UOO 


a op;o 

u.yoz 


/I O 

4.y 


4.0 


A A 
4U 


A OOOOA7 

u.yyyyu ( 


1 aqc; 

iyyo 





T3n n/itr Q TTT 

Jougey-o-iii 


6t ■ TO 
1j1 - Lb 


U.000 


a qos 
U.ozo 


A A7S 

U.U / 


U.UOO 


A G7ft. 
U.O / 


1/11 
r4. 1 


A Q 
4.0 


yo 


A 000/170 

u.yyy4 1 y 


Iyyo 


ft 



Lroesgen-1 


3 TT„ 1 T C 

rie + Lb 


a con 
U.OzU 


0.2(4 


A A7/1 
U.U(4 


A A A O 

U.U42 


a nftft 

U.yoo 


D.O 


C A 
O.U 


00 


A AAAA1 ft 

u.yyyyio 


1 AOft 

lyoo 


7 
f 


Goesgen-II 


3tt„ 1 T c 

rie + .Lb 


a p;c/i 
U.004 


n 00s 
U.zyo 


A OAS 

U.UOO 


a n^A 
U.UOU 


A OOO 

u.yyz 


D.O 


A 

o.U 


A K 
40 


A OOOSSQ 

u.yyyooo 


iyoo 


o 
O 


Goesgen-III 


3tt_, 1 T c 

rie + Lb 


U.04o 


A QOA 

u.ozy 


A A7A 
U.U 1 U 


U.Uoo 


A HOC 

u.yzo 


/ .0 


ft A 

O.U 


00 


A AAA7Cft 

u. yyy / 00 


1 AQft 

iyoo 


Q 

y 


TT T 


ne t" Ijo 


~l 








n ano 

U.OUZ 


q 

y.o 


O.U 


Q 

y 


a qqqqqp: 

u.yyyyyo 


1 QS1 

iyoi 


1 a 
1U 


Krasn. I 


3tt^ 1 ~DT? 

He + rhj 


~1 








u.yoo 


O.O 


A O 

4.y 


QQ 
OO 


A AAAAQ7 
ij.iJiJiJiJO ( 


1 AQ7 


1 1 
±1 


Krasn. II 


ne -j- t Jjj 


~1 
~1 








u.yoo 


on 9 
zu.o 


/I Q 

4.y 


QO 

yz 


A QAQCI 1 

u.yyyoii 


1 Qft7 

iyo < 


1 


Krasn. Ill 


ne -f- t j_j 


~1 








n Q/17 
u.y4 / 


4.y 


A Q 

4.y 


^7 
/ 


A QQQQ1 O 

u.yyyoiz 


1 QS7 

iyo ( 


lo 


CRT) T 

oivr-i 


C*A T Q 
LiU-Lb 


~1 
~1 








a qe;o 
u.yoz 


0. / 


1 


1 Q 
10 


A OOOOQ1 

u.yyyyoi 


1 AQft. 

iyyo 




ODD TT 
O 1 \ 1 - 1 1 


rM t Q 

LjU-IjO 


~1 
~1 








i.uio 


0.0 


7 


O/l 
Z4 


A QQQQft7 

u.yyyyo * 


iyyo 


15 


ROVN088-1I 


3 He + PE 


0.607 


0.277 


0.074 


0.042 


0.917 


6.9 


5.7 


18 


0.999981 


1988 


16 


ROVN088-2I 


3 He + PE 


0.603 


0.276 


0.076 


0.045 


0.948 


6.9 


5.7 


18 


0.999981 


1988 


17 


ROVN088-1S 


Gd-LS 


0.606 


0.277 


0.074 


0.043 


0.972 


7.8 


7.2 


18 


0.999981 


1988 


18 


ROVN088-2S 


Gd-LS 


0.557 


0.313 


0.076 


0.054 


0.959 


7.8 


7.2 


25 


0.999964 


1988 


19 


ROVN088-3S 


Gd-LS 


0.606 


0.274 


0.074 


0.046 


0.938 


7.2 


7.2 


18 


0.999981 


1988 


20 


Palo Verde 


Gd-LS 


0.60 


0.27 


0.07 


0.06 


0.975 


6.0 


2.7 


750-890 


0.967 


2001 


21 


Chooz 


Gd-LS 


0.496 


0.351 


0.087 


0.066 


0.961 


4.2 


2.7 


998-1115 


0.954 


1999 



TABLE I. Tabulated results of all 23 experiments. Experiments are categorized into different groups with horizontal lines. 
Within each group, the a corr represent the correlated uncertainties among different experiments. This table is an extension of 
Table. II of Ref. H. There are additional correlated uncertainties, since Double Chooz results were anchored to the Bugey-4. 
See the text for more explanations. 
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